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The possibility of detecting the S-wave of the decays r — > Stt + j/t- in the threshold region is explored, with 
emphasis on the sensitivity to the size of the quark antiquark condensate (qq) . 
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1. INTRODUCTION 

Despite a great deal of efforts on both 
the theoretical and phenomenological side (see 
Ref. |l|] and references therein), the mechanism 
of the spontaneous breakdown of chiral symmetry 
(SBxS) in QCD remains still unknown, in partic- 
ular for what concerns the size of the simplest 
order parameter of this breakdown, the quark- 
antiquark condensate (qq) . If most of (quenched) 
lattice simulations support the standard hypoth- 
esis that SBxS is driven by the formation of a 
large (qq) condensate, there are indications that 
dynamical massless flavors could strongly modify 
the situation |^. On one hand some unquenched 
simulations seem to indicate that chiral phase 
transitions could happen for a number of mass- 
less flavors Nf as small as Nf = 40]. On the 
other hand, a recent sum-rule evaluation of the 
Zweig rule violating low energy constant Lgj sug- 
gests that the quark condensate in the SU(3) chi- 
ral limit could be considerably smaller than the 
same quantity in the SU(2) chiral limit W. On 
the experimental side, on-going (E865 at BNL, 
DIRAC at CERN) and forthcoming (KLOE at 
Dafjjne) facilities, will allow to pin down (qq) in 
the SU(2) chiral limit from low energy tttt observ- 
ables p. Due to the enormous r-data sample 
already accumulated (~ lO"^ t pairs at CLEO) 
and to the improvements expected in the near 
future § (CLEOIII-CESR, BaBar,. . . ), it is nat- 
ural to explore the possibility of measuring (qq) 
from T decays, which would provide an indepen- 
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dent cross-check of determinations from tttt scat- 
tering. 

The sensitivity to (qq) in the exclusive decays 
r ^ Stt + i^^ is contained in the S-wave. The 
latter, being proportional to the light quark mass 
m = (to„ -|- TOd)/2, is small compared to the P- 
wave contribution, which, moreover, is enhanced 
by the resonance ai. However it turns out that, 
close to threshold, the P-wave is kinematically 
suppressed relatively to the S-wave, thus allowing 
the latter to appear through sizeable azimuthal 
left-right asymmetries. 

2. KINEMATICS 

There are two different charge modes in the 
decay r — > 3tt + Vt, the 27r"7r~ mode and the 
all charged mode, 27r~7r+, whose hadronic matrix 
elements, H"°-{pi,p2,P3) and H--+{pi,p2,P3), 



ijO"- = (^"(pi)vr°(P2)^-(P3)|A;|0), 
H--+ = (^-(pi)7r-(p2)^+(P3)|A7|0), 



(1) 



with At^ — U7^75(i. In the isospin limit, mu = 
nid = TTi, the most general Lorentz decomposi- 
tion of each matrix element involves three form 
factors. 



<'-^lX''bl,P2,P3) 



+ V2^F^'%Pl,P2,P3) + Vi^FriPl,P2,P3),m 

where hfs (hadronic final state) stands for 00— or 
h and 



vr^Pt~P^3-^^^^Q^ 

^2 ~ P2 Pa Q2 V I 

vi'^pt+p'2+p'i^Q^- 



(3) 



For both the hadronic final states, Bose symme- 
try requires that ^2(^1,^2, Pa) = -Fi(P2,Pi,P3)- 
The form factors Fi {F2) and F4 correspond re- 
spectively to J = 1 and J — 0, where J is the 
total angular momentum of the hadronic system 
in the hadronic rest frame. The differential decay 
rate is given by 



dT{r ^ iTTi^r) = ^^\M\^d<i>,, 



(4) 



where A4 is the matrix element of the electroweak 
interaction and d4>4 is the invariant phase space 
of four particles. As shown explicitly by Kuhn 
and Mirkes W , the matrix element squared can be 
expressed in terms of 9 independent leptonic and 
hadronic structure functions Lx and Wx- The 
hadronic structure functions Wx only depend on 
the hadronic variables, Q^ and the Dalitz plot 
variables si = (p2 + PsY and S2 = {pi +^3)^, 
while all the angular dependence is relegated in 
the leptonic structure functions Lx ■ Four of the 
hadronic structure functions correspond to the 
square of the spin 1 part of the hadronic matrix 
element, 

Wa = {xl+xl)\Fif + {xl+xl)\F2f 

-^2(2:12:2 - xl)Re{FiF;), 
Wc = {xl~xl)\F,\^ + {xl-xl)\F2\' 

+2{xiX2 + xl)Re{FiF;), 

Wd = 2[xiX3\Fi\^-X2X3\F2\'' 

+X3{x2 ~ xi)Re{FiF;)] , 
We = -2a;3(a;i+a:2)Im(i^iF2*), 

one is the square of the spin component 



(5) 



WsA = Q^\Fi\^. (6) 

and the remaining ones are the interference be- 
tween the spin and spin 1 components 



WsB = 2^xiRe{FiFi) + (pi ^ P2), 
Wsc = ~2,/Q^xilm{FiFi) + {pi ^ P2), 
WsD = 2^x3ReiFiFi) - {pi ^^2), 

WsE = -2y^X3lTniFiFl) - (pi ^ P2), 



(7) 



frameg, xi = Vf, X2 = V2 and X3 = V^: they 
vanish at the production threshold. Their explicit 
expressions in terms of Lorentz invariant quanti- 
ties can be found in Ref . [0 . As it is clear from the 
explicit expressions of the Wx , the purely spin 1 
structure functions are suppressed, at threshold, 
by two powers of the Xi's, the interferences of 
spin 1 and spin by only one power whereas the 
purely spin structure function Wsa is not sup- 
pressed at all. This is the reason for expecting a 
large sensitivity to the S-wave, and then to {qq) 
in the threshold region. 

3. GENERALIZED CHIRAL PERTUR- 
BATION THEORY 

The appropriate tool for studying the depen- 
dence of the form factors on {qq) at low energy is 
the generalized version of chiral perturbation the- 
ory (GxPT) |§. The results of the exphcit cal- 
culation of the form factors at one-loop level can 
be found in Ref. Q. Taken in the standard yPT 
limit iQ , they reproduce the results of Ref. Q . 
It is convenient to parametrize the dependence 
on {qq) of the form factors of t -^ 3t: + i>t in 
terms of two parameters a and /3, related to 
the elastic tttt scattering amplitude. At leading 
ordeijj they correspond respectively to the am- 
plitude and the slope at the symmetrical point 
s = t = u = 4/3A/2, 






1-^ 



(8) 



In fact a and /? are correlated quantities: they 
must satisfy the so-called Morgan-Shaw univer- 
sal curve, which results from the analysis of Roy 
equations (see Ref. [|2| and references therein). 
All the dependence on {qq) is contained in the 
parameter a: at leading order, it varies between 
1 and 4 if the quark condensate decreases from 
its standard value down to zero, while (3 stays 
always close to 1. The relationship between the 
quark condensate in the SU(2) chiral limit and a 



The Xi are kinematical functions which arc linear 
in the pion three-momenta in the hadronic rest 



^ This system is oriented in such a way that the x-axis is 

along the direction of pa and all the pions fly in the x-y 

plane. 

^ For the definition of a and /3 up to two-loop level see 

Ref. Ili. 




Figure 1. The Gell-Mann-Oakes-Renner ratio 
as function oi a + 2/3. 



and f3 at one-loop level, with the associated theo- 
retical uncertainties, is shown in Fig. n^, where we 
plot the Gell-Mann-Oakes-Renner ratio, 



2m 



X2 



F^M^ rh^O,m,^0 



I™ ,Mi)' 



(9) 



as function of a -f 2/5. Also shown in the figure, 
between the dashed lines, is the experimentally 
allowed range for the combination a + 2/3, as ex- 
tracted from a fit of the two-loop tttt amplitude 
to the Ke4 Rosselet dat a [|l3| . The latter fit yields 
for a and /3 the result |l2|, 

a'^'^P = 2.16 ± 0.86, /3'=''p = 1.074 ± 0.050, (10) 

while the predictions of standard xPT read ||lj] 

a'^' = 1.07 ± 0.01, /3"' = 1.105 ± 0.015. (11) 

4. AZIMUTHAL ASYMMETRIES 

The hadronic structure functions which we are 
interested in are the interferences between the 
spin and spin 1 amplitudes. These structure 
functions contain the dependence on (qq) through 
the S-wave component and, at the same time, 
they are not as small as the purely spin struc- 
ture function, which is suppressed by two powers 
of the light quark mass rh. If the r rest frame can 
not be reconstructed (which is always the case ex- 
cept for r-charm factories), only two of them are 



measurable, Wsb and Wsd- They are related to 
the azimuthal asymmetries, obtained by integrat- 
ing the decay rate over all angles except the az- 
imuthal angle 7 (for the definitions of all angular 
variables see Ref. 0), 



dP = 



Gp 



K 



128Af^ (27r)5 



M2 



n 2 



A/2 + 2g2 



3M2 



xf{Q'n)w{Q')dQ'^, 

Zn 



with W{Q^) defined by 



W{Q^) ^WA + 



MC 



;WSA, 



(12) 



(13) 



M2 + 2Q2 

and the azimuthal distribution, normalized to 1, 

f{Q\l) - l + A2(C[RCOs27-|-C(jDsin27) 
+Ai (Clr cos 7 + Cud sin 7) . (14) 

We have denoted by lowercase letters wx the 
corresponding structure functions Wx integrated 
over the whole Dalitz plot. The asymmetry coef- 
ficients Clj^, C(ji3, Clr and Cud in Eq. ( [l4|) are 
related to the Kuhn and Mirkes' structure func- 
tions by the relationgQ] 



Clr 

Cud 
Clr 
Cud 



1 _ w_ 

^ Mi 



3Mi 



wc_ 



1 



M2+2Q2 W 

TiJ ) TTfpiQ^ w ■ 

_ 7r 3A'/^ WSB 

4 Mj+2Q^ W ' 
7r 3M^ WSD 
4 Mj+2Q^ W ■ 



(15) 



The coefficients Ai and A2 J^ are kinematical 
functions resulting from the integration over the 
T-decay angle. They depend on Q^ and on the r 
velocity. In the limit of ultrarelativistic r (rele- 
vant e.g. for CLEO), they take the form 



Ai(Q') = 



M^ -Q^ + 2M2g2 log -S 



M? 



(Af2 ^ Q2y 



A2(Q2) = -2 + 3 



A/2 + Q2 
A/2 - Q2 



Ai(Q2). 



(16) 



The explicit calculation of the form factors at one- 
loop level shows that the structure function wsd 

* We are neglecting the polarization of the r's, which is 
justified provided that the Z exchange can be neglected 
(far from the Z peak). 
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Figure 2. The integrated left right azimuthal 
asymraetry for the two charge modes. 



vanishes at leading order, unlike wsb- The latter 
governs the left-right asymmetry. In Fig. H we 
plot the quantity AN{Q^) = \N l{Q^) - N r{Q'^)\ 
for a total number of lO'' r-pairs produced, where 
N{Q'^) denotes the total number of events, inte- 
grated from threshold up to Q^, and the subscript 
L (R) refers to events with tt/2 < 7 < 37r/2 (or 
the complementary interval). In order to study 
the sensitivity to (qq) this quantity is plotted for 
three values of a and /3: for each charge mode 
the lower band corresponds to the standard pre- 
dictions (O) , the middle band to the central "ex- 
perimental" values (pJ]| ) and the upper band to 
a = 4, P = 1.16. For both the charge modes 
AA^(0.35GeV^) changes by a factors 3 if a varies 
between its standard value and 4. Quantitatively 
the effect is larger for the all charged mode, large 
enough to be eventually detected with a sample 
of 10^ T-pairs. Each band represents the theo- 
retical uncertainties coming from the low energy 
constants, but not from higher order chiral correc- 
tions. The influence of the latters is estimated by 
studying the convergence of the chiral series. The 



lines immediately below the bands correspond to 
the result at tree level. We see that the one- 
loop corrections amount approximately to 30% 
at Q^ = 0.4 GeV^. Therefore, supposing a ge- 
ometrical behavior of the xPT series, we expect 
the two-loop corrections to modify the results by 
no more that 10%. 

5. CONCLUSIONS 

We have identified the decays r -^ 37r -(- z^t- in 
the threshold region as a promising observable 
for measuring the strength of the quark conden- 
sation in QCD. A total statistics of 10^ r-pairs 
seems sufficient to provide an interesting cross- 
check of future determinations of (qq) from tttt 
observables. 
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